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INERTIAL FORMS OF NAVIER-STOKES 
EQUATIONS ON THE SPHERE 

Roger Temam^'^ & Shouhong Wang^ 



Abstract. The Navier— Stokes equations for incompressible flows past a two-dimensional 
sphere are considered in this article. The existence of an inertial form of the equations is 
established. Furthermore for the first time for fluid equations, we derive an upper bound 
on the dimension of the differential system (inertial manifold) which fully reproduces the 
infinite dimensional dynamics. This bound is expressed in terms of Grashof Numbers. 



Introduction 

Inertial manifolds are smooth finite dimensional manifolds that attract all orbits of 
a dissipative dynamical system at an exponential rate. By restricting the dynamical 
system to this manifold we obtain a finite dimensional one called the inertial form of 
the system. The inertial form of a system has exactly the same dynamics as the initial 
one, due in particular to the so-called asymptotic completeness property (cf. [CENT], 
[FSTi]). 

The introduction of inertial manifolds in [FST] and their name were motivated by 
the hope that inertial manifolds would exist for the Navier-Stokes equations (at least in 
dimension two) and that they would be able to describe the inertial regimes of turbulence. 
This project has been hampered by a restrictive condition for the existence of inertial 
manifolds which appeared in [FST] and in all subsequent works on inertial manifolds, 
namely the spectral gap condition. This condition is not satisfied by the Navier-Stokes 
equations, even in space dimension two. A breakthrough was recently made by M. Kwak 
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[Kl] who embedded the Navier-Stokes equations on a torus (space periodic flows) into a 
reaction diffusion system for which the spectral gap condition is satisfied when the ratio 
of the two periods is a rational number. In this case Kwak obtains, in this manner, an 
inertial form for the Navier-Stokes equations. 

In this article we want to address the case, very important for meteorology, of the 
Navier-Stokes equations past the sphere. We use the same methods as in [Kl] but our 
treatment is totally different because, in particular, we encounter here the difficulties 
related to the geometry. We are able as in [Kl] to reduce these equations to a reaction- 
diffusion system for which the spectral gap condition is satisfied. Furthermore we are 
able in our case to derive an upper bound on the dimension of the inertial manifold in 
terms of the Grashof number. Although this bound is high, it is still polynomial (and 
not exponential) with respect to the Grashof number. 

The article is organized as follows. Section 1 describes the Navier-Stokes equations on 
5"^. In Section 2 we introduce the embedded system and study its relation to the initial 
one. Absorbing sets and attractors for the embedded system are studied in Section 3. 
Inertial manifolds are obtained in Section 4. The dimension of the inertial manifold 
is estimated in Section 5. Finally we conclude in Section 6 with some complementary 
remarks on fiows past a torus and on the barotropic equation of the atmosphere. 



Consider the spherical coordinate system (x^, x^) = {9, ip) on (S^, which 6 G (0, tt) the 
colatitude and (/? e (0, 27r) the longitude. The Riemannian metric is given by 



with its inverse denoted by {g^^)- For geometric notations, see e.g. [A] and for, Navier- 
Stokes flows. Sec. 3.4, Chapt. Ill in [Tl]. 

For simplicity, we use u = u^-^ = {u^) to denote a vector field on S'^. We can then 
write the 2D Navier-Stokes equations on S"^ as follows: 



where the covariant derivative V^iU, the divergence div, the gradient V and the Laplace- 
Beltrami operator A are defined by 



1. Navier-Stokes Equations on S'^ 




(1.1) 




(1.2) 



< 
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(1.4) 



As we discussed in the introduction, we consider in this article the stream function 
formulation of the Navier-Stokes equations. To this end, we define the curl operators for 
both scalar and vector functions as follows: 

{curl i/j = —n X V'i/', 
curl u = — div (n x it), 

where n is the outward normal vector, i/; is a scalar function, and w is a vector function 
on S^. Then we can show that 

jv.„^v(f 

I Au = V(divtt) + {u X n) curl u. 

Now we set 

u = curl t/j, 
( = curl u = — Ai/j, 

where ip is the stream function. For simplicity, let 

(1.5) Ao = -A:Ho^ Ho, 

be the minus Laplace-Beltrami operator on S"^ for scalar functions with inverse Aq^ : 
Ho — > Ho- The function space Ho here and the function space Vq used later on are 
defined by 

Ho^L\S')/R^{CeL\S')\[ 
Vo = H\S\ 

where H'{S'^) = H^S"^) /R for any s G M. Then (1.4) can be written as 

(1.7) C = Ao^, V = ^0 'C, u = curl {A'^\). 

We then infer from (1.1) the following stream function formulation of the 2D Navier- 
Stokes equations on the sphere: 

( Ct + div [Cu] + lyAoC = F = curl /, 

(1.8) }u = cm\{A^\), 

[ C\t=o = Co- 

Obviously, the operator is an unbounded, self-adjoint linear operator with compact 
inverse, and with domain D{Ao) = H'^{S'^). Moreover, using the spectrum of we can 
define its power ^ ^ ^ ^^d 

A'o : H'+^ H', Vs e R 

is an isomorphism. 



(1.6) 
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2. The Embedded System 



(2.2) 



2.1. The Equations for VC and (u. 

Let Co e D{Ao) and ({t) be a solution of (1.5). We define an injection map J : Vb — > 

by 

(2.1) J{C) = {C,v,w) 

where 

'v = WC, 

t(; = Ccurl {Aq\) = Cu. 

The Hilbert space H used here and the Hilbert space V which will be used later, are 
defined by (cf. Remark 2.1 below) 

( H = HoxHixH2, 
H^ = VH\S^),H2 = L\TS^), 
V = VoxViX V2, 
Vi = VH^{S^), V2 = H\TS^). 

According to Hodge decomposition theorem (cf. [A]), we have 

C°°{TS^) = {V0 I e C°°{S^)} e {curlV' I V e C°°{S^)}, 
which implies that 

(2.3) H\TS'^) = VH'+\S'^) © cmlH'+\S'^), Vs > 0. 

Therefore, Hi and Vi in (2.2) are wcU-dcfincd. 

Moreover, for simplicity, we need to define two linear unbounded operators Ai on Hi 
{i = 1,2) as follows: 

Ai = -A:Hi^Hi, 
A2 = -A:H2^ H2, 

the operator Ai being the restriction of —A in Hi = 'S7H^{S'^) with domain D(Ai) = 
\7H'^{S'^), and the operator A2 being the minus Laplacian —A in H2 = L'^{TS'^) with 
domain ^(^2) = H^{TS^). 
Then we have 

Lemma 2.1. The functions J{() defined by (2.1) satisfies the following system of 
equations 

Ct+ diyw + iyAoC = F, 

(2.4) {vt + ^ (divw) + vAiv = VF, 
wt + i'A2W = Fu + Cf - 2vTr[v ® Vu] - {v • u)u - ( curl [Aq-^{v ■ u)], 
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where u = curl {Aq^Q. 

Proof. Obviously, (2.4)i is the same as (1.8). Applying the gradient operator to 

(2.4)1, we obtain easily (2.4)2. 

By definition, and using (1.8) and ut — curl {Aq^C^i), we find 

wt = Ctu + Cut 

= {F-diy[Cu]-iyAoC}u + Ccm\{A-^[F-diy[Cu]-iyAoC]}. 
We then infer from the identity 

A{av) = {Aa)v + aAv + 2Tr((Va) <S> Vv), 

that 

Aw = {AC)u + (Au + 2Tr((AC) Vu). 

Therefore, 

Wt — vAw 

= Fu + Cf -udiy iCu) - C curl [^-^(div {(u))] - 2iyTr[{WC) ® (V«)], 

which implies (2.4)3, using div {(u) = (divu + (VC) ■ u = v ■ u. The proof is complete. 
□ 

2.2. The Embeded System. 

Equations (2.4) by themselves do not seem to be dissipative. Therefore in order to 
capture the dynamics of the NSE (1.1) or (1.5), we shall modify (2.3) in such a way that 
the resulting system satisfies the following properties: 

(a) The dynamics on the section J(Vb) is unchanged, 

(b) The modified system is dissipative, and 

(c) All the solutions of the modified system will converge to the section J(Vo) as 
t ^ 00. 

More previsely, we propose the following modified system of equations: 

' Ct + J^AoC + dw w = F, 
vt + lyAiv + V(div w) = VF, 
^^'"^"^ Wt + VA2W — Fu — Cf + i'V ' '^)'^ + C curl [Aq ""^(f ■ u)] 

+ 2vTr[v ® Vu] + kv[l + J^'^HCHhV^]!^ - iD) = 0, 

where w =^ C, curl [Aq^Q. The system (2.5) is obtained from (2.4) by adding the 
underlined term; the positive constant k will be specified later on. Generally speaking, 
5 7^ V, w ^ w{v = VC, w as before). 
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For simplicity, we write Z = {(, v, w). Then we have the following theorems about the 
existence and properties of solutions of the embedded system (2.5). 

Theorem 2.1. We assume that F is given in Hq, and Zq = {Cvo^wq) is given in 
Vq X Hi X H2. Then there exists a unique solution Z = Z(t) of (2.5) with initial value 
Z(0) = Zq, such that 

(2.6) Z e l2(0, T; D{Aq) xViX V2) n C([0, T]; Vb x i^i x i^a), V T > 0. 
Moreover, if f e Vq and Zq e D{Aq) x Vi x V2, then 




Z e L2(0, T; D{AI/^) X D{Ai) X D{A2)) D C([0, T]; L>(Ao) x Vi x V2), 
ZteL^{0,T;Vi xHi x H2). 



Theorem 2.2. Let ({t) be a solution of (1.8) with initial value C(0) = Co £ D{Ao), 

and let f E Vq. Then Z{t) = J{({t)) is a solution of (2.5) with initial value Z{0) = J{zq) 
for all t > 0. Conversely, if Z{t) = {C{t),v{t),w{t)) is a solution of (2.5) with Z{0) = 
J(Co) e D{Ao) xViX V2, then C(t) is a solution of (1.8) fort>0 and Z{t) = J{C{t)). 
In particular, J{D{Aq) x Vi x V2) is a positively invariant set for (2.5). 

Theorem 2.3. For any Zq = {Co,vq,wq) e Vq x Hi x H2, 

(2.8) \v{t) - v{t)f + u-^\w{t) - w{t)f < e'^WyQ - vq^ + u'^lwo - wof), \/ t > 0. 

Consequently, 

lim |Z(t)- J(C(t))| = 0. 

t — >oo 

The proof of Theorem 2.2 is easy. We only have to prove Theorems 2.1 and 2.3. 

Proof of Theorem 2.1. First of all, we solve the equations (2.5) 1 and (2.9)-(2.10) 
below for the unknown functions (C, v — v,w — w), by Galcrkin method. Then we obtain 
the solution Z = {(,v,w) of (2.5). Since the Galerkin method is a standard procedure, 
here we only present some formal a priori estimates, which will also be used later for the 
existence of attractors of the embedded system (2.5). 

Step 1. Equations for — v and w — w. We infer from (2.5) 1 that 

vt + vAiv + V(div w) = VF, 

which implies that 



(2.9) 



{v — v)t + vAi {v — v) = 0. 
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We also have, as in the proof of Lemma 2.1 and since div w = v ■ u, 

wt + i'A2W — Fu — Cf ~ 2vTr[v ® Vtt] — (div w)u — C, curl [Aq ""^(div w)] 

Therefore 

(2.10) {w-w)t + -i^A2{w-w) 

-\- {v ■ u — div w)u + C curl [Aq ""^(i; ■ u — div w)] 

+ 2vTr[{v -v)® Vu] + kv[l + i^~^\\C\\m/2\{w - w) = 0. 



Step 2. Energy estimates I. We infer from (2.5) i that 

l_d 
2dt' 



i^lCr + H4/'Cr = (i^,C)-(div«;,C) 



< |F|.|C| + |(div («;-«;), 01 + |(div«;, 01 

< |F|-ici + h-*i-iA;/\i 



I.e. 



(2.11) ||c|2 + .|A;/^ci^<^(|i^r + k-^i^). 

By (2.9), we have 

(2.12) ^\v-v\'' + 2v\Ay'^{v-v)\'' = {). 
Moreover, we deduce from (2.10) that 

(2.13) ^^\w - + ua\Al^\tv - w)\^ + kua[l + iy~^\\C\\U/^] ■ \w - 

1 /2 

= —a{{v ■ u — div w)u + C curl [Aq [v ■ u — div w)] 
+ 2vTr[{v -v)® Vu],w - w), 

where we introduce for convenience a constant a = > 0. 
The 1** term in the right hand-side of (2.13) is equal to 

—a{[{v — v) ■ u + div (w — w)]u, w — w) 

< ca\\u\\'j^8 ■ \w — w\ ■ \Ay^{v — v)\ + ca\\u\\L'^ ■ \w — w\ ■ \A2^'^{w — w)\ 

< (since llClk^ - IIAV'lk^ - ~ HVV'IIh^+i) 

< ^{a'\C\' + a||C||?,v.) ■\w-w\' + l\Ay\v - v)\' + ''-^\A}I\w - w) 



8 R. TEMAM & S. WANG 

The 2"*^ term in the right hand side of (2.13) is equal to 

— Q!(Ccurl {^^""^[(v — v) ■ u + div {w — w)]}, w — w) 

< ca[\\C\\l^;My\v - v)\ + \\C\MAy\w -w)\]-\w-w\. 

The 3^'^ term in the right hand side of (2.13) is majorized by 

1/2 ~ 

cz^q;|w — w| • ||C||j:fi/2 {v — v)\. 

Therefore 

(2.14) a^\w - tup + lyalAl^'^iw - w)\'^ + 2ki^-^[l + iy-^\\C\\%^/2] ■ \w - w\'^ 



2 



< ciy-\l + iy-^\\C\\%./2) ■ \w - + iy\AY\v - v) 
Then, with the constant a = v~'^ still at our disposal, we have 

j^WCf + |v - 51" + a\w - wf\ + iv[|^J/'cl' + l^;^'(^ + Oi\A)l''{w - w) 



+ 2kv ^W^y- iMi^i/2j 



Hence for k (absolute constant) large enough, we obtain 

(2.15) + 1^ _ + _ ^|2] ^ v\\A]l\f + \A^I\v - v)f 

+ a\A}l\w - w)f\ + kv-\\ + v-^\\C\\\ir,.\ ■\w-w\ 

^ 2 

V 



<^|F| 



Furthermore we consider the inner product between (2.5) i and ^oC) ^^^d we obtain 
1 d 



2dt 



\Al/\f + u\AoCf = {F, AoO - (div w, AoQ 

< \F\ ■ \AoC\ + I (div {w-w),AoC)\ + I (div w,AoC)\ 
<^\F\^ + + ^l^oCP + Idiv *| ■ |AoC|. 

Since 

Idiv «;| • l^oCl < ^|CM4^'CP + ^l^oCP, 
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we obtain 



(2.16) + H^oCI^ <"-\Ff+ '-^\Al/\w - w)\^ + ^ICr • l^/'Cr. 

Finally the combination of (2.12), (2.14) and (2.16) shows that 



d 
di 



- — + \v — f 1^ + a\w — w 



r,|2 



+ U 



1 

2ci 



— \AoC\'+\Al/\v-v)\' + ^\Al/\w-w)\' 



+ 2ku-'[l + i^-^C\\W]-\w- 



w\ 



-4||Al|4 



).\w-w\' + ^\C\'-\Al/\\^ + ^\Ff, 



I.e. 



(2.17) 



d 
di 



a\w — w\ 



w — w)\'^ 



+ /c^.-Mi + ^-'IKII^i/2]-k-«;r 
< -s\C\' ■ lAl^W' + -\F\\ 

By Gronwall's inequality, we obtain from (2.15) and (2.17) that 

(2.18) {C,v-v,w-w) e L^{0,T;D{Ao)xVixV2)nL°°{0,T;VoxHixH2), V T > 0. 
Then from equation (2.5)i and (2.9)-(2.10), we can obtain a priori estimates for 

{C,v-v,w-w)t e L2(0,T;i?o x V( x V^), V T > 0. 

Hence 

(2.19) {C,v-v,w-w) e C([0,T];Vb x ifi x H2), V T > 0. 
On the other hand, by definition, we can prove that 

' \v\ = IVCI < lA'J^l, 
\aY^v\ < \AoC\, 
\w\ = IC curl (^-^C)! < 
[\Al/'w\<\AoC\-\Al/\\. 



(2.20) 
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Then the combination of (2.18)-(2.20) proves the existence of a solution Z = (ry, f , w) of 
(2.5) satisfying (2.6). 

Step 3. Energy estimates II. From (2.5) i, we also have 

^^l^oCP + H^^CP = {F,AlO - (div w,AlO 

< \aI/'f\ . \aI/\\ + |(div {w - w), AlO\ + |(div w, AlO\ 

< {\aI/'f\ + |V(div (w - w))\ + |V(VC ■ curl {A-\m ■ |Af CI 

< ^\Ai/'F\^ + ^\A,{w - w)\^ + ^waw ■ + ^l^f CP, 

i.e. 

(2.21) ||AoCP + < ^l^y^Fp + '-^\A,iw - w)\' + ^llClP^v. • \M\'- 
By (2.9), we have 

(2.22) ;^l^y'(^ - ^)P + MAiiv -v)\^ = 0. 
Moreover from (2.10) we obtain that 

(2.23) 

%j^A}l\w - w)\^ + va\A^{w - w)\^ + kv'^X + ^-^||C|ll,v.] • \A'-J\^ - 

= —(xiiy ■ u — div w)u + ( curl [A(^^(t; ■ u — div w)] 
+ 2vTr[{v -v)® V-u], ^2(10 - w)). 

the 1** term in the right hand side of (2.23) is equal to 

— a(\{y — v) • u + div {w — w)]u, A2{w — w)) 

< c«{|Cp • \A\'\v - v)\ + IICII^v. • \Ay\w -w)\]- \A2{w - w)\. 

The 2'^'^ term in the right hand side of (2.23) is equal to 

— a{C, curl {Aq ""^[(t; — v) ■ u -\- div {w — w)]}, A2{w — w)) 

< ca[||C||^^i/2|C| • \A\'\v - v)\ + WCWmfA^ - w\] ■ \A2{w - w)\. 

The 3'^'^ term in the right hand side of (2.23) is majorized by 

cv-^\A2{w-w)\.\\C\\Hy.\Ay\v-v)\. 
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Therefore 



(2.24) 



w 



ca , 



r.M2 



|2 I 4I/2, 



Finally (2.21)-(2.22) and (2.24) imply that 
(2.25) 



d 
dt 



-^lAoCl' + \Ay\v -v)\^ + a\Al/\w - w 
ZC2 



r.M2 



+ 2ku-^[l + u-^aW]-\Al^\w-w)\' 



-L\AoCf + \Ay'{v-v)\' + a\Al/\ 
ZC2 



w — w] 



Using Gronwall's inequality as in Step 2, we can prove that if Zq G D{Ao) x Vi x V2, 
then Z satisfies (2.7). 

The proof of the uniqueness of solution is easy, and is omitted. 

Theorem 2.1 is proved. □ 



Proof of Theorem 2.3. We infer from (2.12) and (2.14) that 

d 

— [|f — + a\w — wp] < — z^[|A]^^^(f — v)'^ + a\A2^'^{w — w 

(ajL 

< —i^llv — + a\w — w\^], 

which proves (2.8). 



t.m2i 



□ 



3. Absorbing Sets and Attractors 



First of all, Theorem 2.1 indicates that we can define the operators 
(3.1) S(t):Zo — >Z(t), Vt>0. 
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These operators enjoy the standard semigroup properties, and they are continuous oper- 
ators in Vq X Hi X H2 and in D{Ao) x Vi x V2. 

We now prove the existence of an absorbing set of the semigroup {T,{t)}t>o in Vb x 
Hi X H2. T othis end, by GronwaU's inequahty, we infer from (2.15) that 

|C(t)|' + |(^ - v){t)\^ + a\{w- < c|F|2 + e-^*(|Co|' + l^^o - ^o|' + «ko - *o|'), 

for any t > 0. Then for any R > 0, if Zq = {CotVq,wo) G BvoxHixH2{^i R)i the baU in 
Vq x i^i X H2 centered at with radius R, there exists to = to{R) > such that for any 
t>to 



(3.2) 



|C(t)P + \{v - v){t)\' + a\{w- wmf < c\F\^ 

[^'uTm? + \A'\v - v){t)\^ + a\Ay\w - w){t)\^]dt < c\F\ 



Now we return to (2.17). By the uniform Gronwall lemma (cf. Lemma 1.1 in [Tl]), 
we obtain that 

(3.3) |A;/'CWI' + |(^ - v){t)f + \{w - w){t)f < c|F|2exp(c|F|4), V t > to + 1. 
By (2.20), (3.3) shows that 

(3.4) \Al/\{t)f + \v{t)f + \w{t)f < c\ff exp(c|F|4) + c\F\^ exp{c\Ff), V t > to + 1- 
Now we define the right hand side of (3.4) as Pq, then we conclude that 

(3.5) msvo 

In other words, By^y^HixH2{^i Po) is an absorbing set in Vq x ffi x H2 for the semigroup 
Moreover, from (2.17), we also have 



(3.6) J [|AoC(t)P + - v){t)\^ + a\Al/\w - wX^dt 

< c|F|2 + cplil + \F\^), V t > to + 1. 



Now we want to establish the existence of an absorbing set in D{Ao) x Vi x V2. We 
apply the uniform Gronwall lemma to (2.25). Then we obtain 

(3.7) \AoC{t)\' + \Al/\v - vm\' + a\Al/\w - w){t)\^ < cid^/^Fl), V t > to + 2, 
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1/2 1/2 

where ci(|Aq' F\) is a function of \Aq F\. On the other hand, it is obvious that 




Ay^v\<c\AoCl 
a1^^w\ < c\AoCf. 



Therefore 

(3.9) \AoC{t)f + \Ay^v{t)f + a\Al/^w{t)f < C2{\aI^^F\), \/t>to + 2. 

1 /2 1 /2 

where C2(|Ao F\) is a function of \Aq F\. This shows that the ball -BD(yio)xVi (0; Pi) 
in D{Ao) X Vi X V2 is an absorbing set in D{Ao) x Vi x V2, and that the operators S(t) are 
uniformly compact for t large. Here pf = C2{\Aq^'^F\). Hence by Theorem 1.1, Chapter 
I of [Tl] we promptly obtain that 

Theorem 3.1. There exists a unique global attractor A C D{Ao) x Vi x V2 of the 
embedded system (2.5), such that A is compact, connected and maximal in Vq x Hi x H2 
and D{Ao) xViX V2. 

The following theorem shows that A is exactly the lifting of the attractor ^0 of 
the original Navier-Stokes equations (1.1) or (1.5), via the injection map J : Vq — > 
Vox Hi X H2 defined by (2.2). 

Theorem 3.2. We have 

(3.10) A=J{Ao). 



4. Inertial Manifolds 



4.1. Preliminaries. 

The embedded system can be rewritten in the following concise form 



(4.1) —+AZ + R{Z) = 0, 

where the linear and nonlinear operators A, R are defined by 

/ lyAoC + diYW \ fRo{Z)\ 

(4.2) AZ= i^Aiv + V{div w) , RZ = [ Ri{Z) \ , 

\ l^A2W J \R2{Z)I 
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where Ro{Z) = -F, Ri{Z) = -VF and 

R2{Z) = -Fu-Cf+{v-u)u+C curl [A-^{v-u)]+2i^Tr[v^Vu]+ki^[l+i^-'^\\C\\%i/2]{w-w). 

Then we can obtain easily that D{A) — D{Ao) x D{Ai) x D{A2), and D{A) is dense 
in H. 

For technical reasons, we need to use the following equivalent norm \ • \h iu H: 

\z\HH\C\' + \v\' + o^H'y^', 

the corresponding inner product being denoted by {■,-)h', here a = Then we have 

the following Garding inequality 

(4.3) {AZ, Z)h > ^{l^rcr + lA^'vf + o^\A^'w\% 

which can be proved easily by taking the following inequality into consideration: 
(V(div w),v) = («;, V(div v)) = {w,Ai,v) < llA^^v]'^ + 1\aI^'^w\'^. 

Moreover, as in [Kl], we can prove that A is a sectorial operator, and A~^ : H ^ H isa 
compact linear operator. Therefore the spectrum a (A) consists of a countable number of 
eigenvalues with no finite accumulation points and each eigenvalue has finite multiplicity. 

Indeed we have 

Lemma 4.1 (1). The operator A has only the eigenvalues un{n + l)(n = 1, 2, . . . ). 

(2). The eigenvalue fnin + 1) has at most multiplicity 8(2n + 1). 

Proof. (1). We can easily prove that if A is an eigenvalue of A, then A has to be an 
eigenvalue of either vAq or vAi, or i'A2. Then it suffices to prove that all the eigenvalues 
of Aq, Ai and A2 are the numbers {n{n + l)|n — 1,2,...} and count multiplicities. 
According to [CH], the eigenvalues of Aq are the numbers {n(n + l)|n = 1,2,...}, and 
the eigenfunctions corresponding to n(n + 1) are the 2n + 1 spherical harmonics, denoted 
Yn,h{h = 0, ±1, . . . , ±n). Thanks to the Hodge decomposition (2.4), we can easily prove 
that \/Yn^h{h = 0, ±1, . . . , ±n) are eigenvectors of Ai corresponding to the eigenvalue 
n{n + 1), and that {VYn^hlh = 0, ±1, . . . , in, n = 1,2,...} is a complete basis of Hi. 
Therefore, the operator Ai has only eigenvalues {n(n + l)|n = 1,2,...}. 

Similarly, we can also prove that the operator A2 has only eigenvalues {n{n + l)|n = 
1,2,...} with corresponding eigenvectors {Vl^n./i + curl Yn^h'\h, h' = 0, ±1, . . . , ±n}. 

(2). Set A = un{n + 1), and 

fAox Bo 
Ax = A-XI=\ A,x B, 

V A2X 
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where Aix = vAi — XI {i = 0, 1, 2), = div and Bi = V(div) and all unwritten elements 
of the matrix vanish. Then by direct computation, we can obtain that 



= Al^ Bf^ 

where Bf^ = ^a^f + BiA^^\i = 0, 1. 

Now we suppose that A'^Z = for Z = {(, v, w). It follows immediately that A2)^w — 
in which case w = or A is an eigenvalue of A2 and w an associated eigenvector. In the 
first case we infer that ^qaC — 0; ^ia^ ~ ^' hence A is one of the numbers vn{n + 1) and 
C and/or v are associated eigenvectors of Aq,Ai. In the second case A is again one of the 
numbers vn{n + l) and w is an associated eigenvector of ^2; w = VYi + curl I25 ^^i and 
Y2 being linear combinations of the Yn^h{h = 0, ±1, . . . , ±n). We infer that 

B^^^w = AoxBow + B0A2XW 

= AoxdiY[VYi + curl Y2] 

= Aoxdiv VYi 

= Aoa(AYi) 

= -n{n + l)AoxYi 

= 0. 

(2) (k) 

Similarly, Bl 'w = and by iteration B^ '^w = 0, for i = 0, 1 and all A; > 2. Hence we 
obtain that 

ker Al = {(C, v, w) \ AqxC = 0, Aixv = 0, A2XW = 0}, 

and this is the generalized eigenspace of A with respect to A. In summary, A has at most 
the multiplicity 2(2n + 1 + 2n + 1 + 2(2n + 1)) = 8(2n + 1). 

The proof is complete. □ 

To prove the existence of inertial manifold for the embedded system, it is common to 
truncate the nonlinear terms R{Z). To this end, we need 

Lemma 4.2. If F e D(Ao), then A is bounded in D(A). 

3 /2 

Proof. It is easy to see that if F e D{Ao), then ^0 is bounded in D{Aq ). Then 
Theorem 3.2 implies that A is bounded in D{A). □ 

Now let p be a positive number such that 

\AZ\h <^,yZeA. 
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We consider the following prepared equations of the embedded system: 

(4.4) _ + AZ + i?^(Z) = 0. 
Here 

(4.5) i?^(Z) = y)p{\AZ\H)R{Z), V Z e D{A), 

where r\p{s) = r\{sj p) and is a C°° function from into [0, 1] such that 

r 1 for < s < 1, 
^(^) = l0 for.>2, 

and sup^>o \r]\s)\ < 2. 

Lemma 4.3. is a globally bounded operator from D{A) into itself: 

(4.6) sup \ARr,{Z)\H < sup \AR{Z)\h = Mi. 

ZeD{A) \AZ\h<2p 

Lemma 4.4. is a globally Lipschitz operator from D{A) ito itself: 

(4.7) \A{R^{Z^) - Rr,{Z2))\H < M2\A{Zi - Z2)\h. 

These lemmas are easy. Explicit values of Mi and M2 are needed and given below. 

4.2. Inertial Manifolds. 

Now we let 

< Ai < A2 < . . . < Atv < Atv+i < . . . 

denote the eigenvalues of A repeated according to their multipliticity. Let P = P/v be 
the projection associated with the eigenvalues {Xn\n < N}, and Q = Qn = I — P- Then 
PH is the direct sum of the eigenspaces corresponding to the eigenvalues A^ < A at, and 
QH is the direct sum of the eigenspaces corresponding to the eigenvalues A„ > Aat+i- It 
is easy to see that 

(4.8) PHcD{A), H^PH®QH, 

and PH, QH are invariant under A. 

The inertial manifold for equation (4.4) is obtained in the form Ai = Graphs, i.e. as 
the graph of a Lipschitz mapping 



(4.9) 



$ : PH ^QHnD{A), 
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such that 

f supp$ C{ye PH, \Ay\H < 2p}, 

(4.10) \ \A^{y)\H<y\ ^yePH, 

[ - A^y2)\H < l\A{yi - ya)!//, Vyi, y2 G PH. 

Theorem 4.1. We can choose N such that there exists an inertial manifold for the 
prepared equation (4:A) of the embedded system (4.1), which is the graph of a Lipschitz 
mapping (4.9) satisfying (4.10). 

Proof. By Lemma 4.1, it is easy to see that there is an integer A'' > 1 such that 



(4.11) 



Aat+i -Xn = 2z/(n + 1) > 2M2^-^. 

Ml 

Aat+i = iy{n + l)(n + 2) > 



Then the result can be obtained by applying the theory of inertial manifolds for nonself- 
adjoint operators to (4.4) (cf. [DT] and [SY]). The estimates (conditions) in (4.11) are 
those of [DT]; b are arbitrarily chosen, e.g. £ = 6=1/2. □ 

Finally we have 

Theorem 4.2 For f e D{Ao), the essential long-time dynamics of the 2D Navier- 
Stokes equations on the sphere S"^ is completely described by the system of ordinary 
differential equations 

(4.12) ^ + Ay = PF{y + ^y)), m PH, 

where ^ is the Lipschitz map (4.9) obtained in Theorem 4--1 satisfying (4.10). 

5. Dimension of the Inertial Manifolds 

In order to estimate the dimension of the inertial manifold for the embedded system 
given by Theorem 4.1, we need to estimate the constants Mi and M2 in (4.11) and the 
relationship between N and n in (4.11). To this end, we have 

Lemma 5.1. Suppose that N is the largest number such that Ajv = vn{n-\- 1), then 

(5.1) iV<8n(n + 2). 

Proof. By counting the multiplicities, we obtain that 

N < 8[(2 • 1 + 1) + • • • + (2 • n + 1)] = 8n(n + 2). 

□ 
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Lemma 5.2. The dimension N of the inertial manifold satisfies 
(5-2) N<l + max^-^, ^J^^^ 

Proof. By Lemma 5.1, we obtain that 

N < 8(n + l)2, 

which imphes that 
Therefore 

[N 

Xn+1 - Aat = 2zy(n + 1) > i^W y . 
Hence if we choose such that 

8M|(1^ 

then (4.11)i is satisfied. Moreover, it is easy to see that if we choose N such that 

8Mi 



N > 



then (4.11)2 is satisfied. 

The proof is complete. □ 

To estimate Mi and M2 we need to estimate p as follows. 
Lemma 5.3. 

(5.3) \AZ\ <p = cu^[GoG4 + Gl% W Z e A, 

where Go and G4 are generalized Grashof numbers^ defined by 



(5.4) Go=^, G4=l'^"^^ 



9 ) 9 



^No length appears in the definition (5.4) of Go and G4 because we assumed that the radius a of the 
sphere is equal to one. For a 7^ 1, we would define Go and G4 by 
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Proof. Step 1. We estimate the bound for Q & Aq m. D{Aq' ). To this end, we have 
from (1.8) that 

i|i<r + Hvim^ = wo < + 

which imphes that 

(5.5) ||C|^ + .|^V\|^<£!^. 

By Gronwall's inequahty, we obtain that 

c\f? 



(5.6) mv< 



for t large enough. Since the attractor is time invariant, by shifting the time t, (5.6) 
holds true for any t > 0. 

Now we derive as in [T3] (see Sec. 4) the higher order norm estimates: 

= {F, AlO - ( VC • u, AlO 

< \AoF\ ■ \AlC\ + |(I^'(VC) ■ u + D{VO ■ Du + i^(VC) ■ D^u, AlQ\ 

< \AlC\ ■ {WfWm + IKIIhv. ■ Id + IICIIhs/. ■ IICIIhv. + IICIlLa/^} 



2' V 



I.e. 



(5.7) ^^\aTQ? + i^m' < + ^ ^ ^-'[Gl + Gl% 
where G3 is the generalized Grashof number defined by 

(5.8) Gs = ^^1^ (< CGj/^cf ^ by interpolation). 
Using Gronwall's inequality, we can show that 

\Al^^C{t)f < ciy^[Gl + Gl% for t large enough. 
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Therefore 

(5.9) \Al/\\'<cu^[Gl + Gl% VCe A. 

Step 2. We estimate the bounds of in D{aI^'^) and D{Aq] 

= {F,AoO-{VC-u,AoC) 
<c\AoC\-{\\f\\m + \\C\\m/^-\C\} 
<c\AoC\-\\f\\m+c\AoC\'/'-\Cf/' 



which imphes that 



j^\Al/\f + u\AoCf<cu' 



c\C\ 



10 



2 

//I 



10 



^^4/3^2/3 , ^10 



It follows that 

Therefore 
(5.10) 



\Al/\{t)f < cv^ fG^/'G^' + Gl^] , for t large. 



Moreover, we also have 



< l^/^Fl . + \{D'C -u + VC-Du, Al/\) 

<I^o^'CI-{||/IIh^ + IICIIhs/-ICI + IICII 

<i^rci-{ii/ii^^ + i^gcr-icr} 



+ 



,11 



Therefore, we obtain that 

d 



dt 



\AoC\' + i^\Al/\\' < cu' Gl/'Gt^' + Gl^ 



INERTIAL FORMS OF NAVIER-STOKES EQUATIONS ON S'^ 



21 



Hence by Gronwall's inequality, we obtain as before that 



(5.11) 



Step 3. Now we estimate the bound for A = J{Aq). To this end, for any Z 
(C, f , w) e A., we have v = VC and w = ( curl {Aq^Q. It follows that 



(5.12) 
(5.13) 



18i 



\A2W\' < cWCWh 

< (by (5.10) -(5.11)) 

<cu''[Gl/^Gt/^ + Gl''] 
Then we can easily obtain using (5.8) that 



lOi 



12i2 



which shows that 
(5.14) 

The proof is complete. 



\AZ\h < cv^[GoGi + G'q% \/ Z eA. 



□ 



Now we need to estimate Mi and M2. 
Lemma 5.4. 

(5.15) Mi<cv^[GlGl + Gl\ 

Proof. For any Z — {(,v,w) G D{A) such that \AZ\h < 2p, we have 
\AR{Z)\l < 2jy2{|^oi?o(z)|2 + |^ii?i(Z)p + a\A2R2iZ)\^} 



+ 2|div i?2(^)P + 2|V div i?2(^)P 
<2^'{||/lll.3 + 11/111,4} + C|^2i?2(^) 

< 2^'{ii/iii.3 + 11/111,4} + c{ii/iii,3iiciii,i + ii/iii,.iicii^. + M'um^ 



1^ 

.2 



2 

H 

2 ^6 



< 2^^{||/||1,3 + 11/111,4} + c{\\f Wls^ + WfWh^ + ^ 

+ ^^(i + ^-V)^(p2 + ^)} 



< (with (5.3)) 
Therefore, by (4.6) we obtain (5.15). 



□ 
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Lemma 5.5. 

(5.16) M2 < cu[GlGl + Gl% 

Proof. As shown in [Tl] (see (2.19) in Chapter VIII), M2 is given by 

(5.17) 



2Mi 

M2 = + C2p, 



where C2p is the local Lipschitz constant of R defined by 

(5.18) \A{R{Z^) - R{Z2))\h < C2p\AiZi - Z2)\h, V e D{A), \AZ,\h <2p,i = 1, 2. 
Obviously, 

2Mi 



(5.19) 



■»60i 



< ciy[G^oGl + 



Now for any Zi G D{A), \AZi\ii < 2p, i — 1,2, we have 

\A{R{Z,) - RiZ2))\H < |div (i?2(^l) - i?2(^2))| 

+ |V div (i?2(^l) - i?2(^2))| + |^2(i?2(^l) - i?2(^2))| 

< (since A2 : D{A2) — > H2 is an isomorphism ) 

<c\A2{R2iZ^)-R2{Z2))\ 



< \A{Z,-Z2)\h 



V 



+ CZv(l + p )\A2{wi — Wi — {W2 — W2))\ 

+ cu-^\A2W,\ + \A2W^\)m\\%,,, - llCsll^i/O 
< \A{Z, - Z2)\n ■ \vG-, + ^ + cv{\ + 1^"^^) 



+ .-(l + .-V)P+^ + ^,, 



Hence 
(5.20) 



ivGa + 4 + ci^(l + ^~ V') + i^-' (1 - V')P +4 + 4 



C2p < 

<cu[GlGl + Gl% 
Then the combination of (5.17) and (5.19)-(5.20) proves (5.16). 



□ 
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We are in position to state the main theorem in this section, which is now a direct 
consequence of Lemmas 5.2 and 5.4-5.5. 

Theorem 5.1. The dimension N of the inertial manifold given by Theorem 4-i is 
bounded as follows: 



Remark 5.1. Let us mention here the result in [EFNT] concerning exponential 
attractors. Exponential attr actors are (nonsmooth) positively invariant sets that attract 
all orbits at an exponential rate. For general dissipative equations, and in particular for 
the 2D Navier-Stokes equations with general boundary conditions, existence of such sets 
was shown in [EFNT]; furthermore the upper estimate of their dimension can be the 
same as that of the attractor, i.e. cGq. For the equation we consider, the dimension of 
the exponential attractor is much smaller than that of the inertial manifold; however it 
is not a smooth manifold. 



6.1. Remarks on the Space Periodic Case. 

We briefly present here a simplified proof of the result of [Kl] on the existence of 
inertial forms for the 2D Navier-Stokes equations with periodic boundary conditions when 
the ratio of the periods is rational. This proof differs from [Kl] in that the embedded 
reaction-diffusion system is much simpler and the analysis is consequently simplified. 

Let X = {x^,x'^) be the coordinate for the 2D torus = R^/ (27rZ)^. We consider the 
2D Navier-Stokes equations on T^: 



(5.21) 



N <c[GlGl + G^o^f. 



□ 



6. Complementary Results 



(6.1) 




where u = u{x,t) = {u^,u^) is the velocity field and p = p{x,t) is the pressure. 
We define as usual the following function spaces (cf. [Tl]): 



Ho = {ue L2(t2)2| / udx^dx"^ = 0, div u = 0, } 



Vb = i^oni^^(T2)^ 

(6.2) \ H, = {Hof,H2 = L\TY, 

V^ = {Vof,V2 = H\TY, 



H = HoxHiX H^, 
V = VoxVixV2. 
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We also consider the orthogonal projectors 
(6.3) 



Pi : L^{TY Hi, 
P : l2(T2)8 ^ H. 



Now we define an injection map J :Vq ^ H hy setting 
(6.4) J{u) — {u,v,w), 

where 



ei(z = 1,2) being the unit vectors in and a;^ directions respectively, with dual basis 
denoted by {e^,e^}. 

Then we can consider the following embedded system 

' Ut = i^Aou + Po{w'J^ei) = f, 
vt + lyAiv + Pi[w%e, ® ej] = V/, 
Wf + UA2W + u ^ PoTriu ® f ) + PoTr{u ®v)®u 
2 

k=i 

+ hull + J^~^||w||^3/4 + i^~^(||'J^||^3/2 + ||'J^||^l/2)](ty - U0U) = 0. 



(6.5) 



It is easy to see that {u,v = v,w — w} satisfy equations (6.5), the underlined terms 
in (6.5) being equal to 0. Considering then the system (6.5), we can prove exactly as we 
did for the NSE on the sphere S'^, that for k large enough, (6.5) is a dissipative system. 
Then all results in [Kl] can be proved. Especially, there is a global attractor A for the 
embedded system (6.5) given by 



(6.6) A = J{Ao), 

Ao being the global attractor of the original Navier-Stokes equations (6.1). Moreover, 
by studying the prepared system of (6.5) as (4.4) one can prove the existence of inertial 
manifolds for the prepared system of (6.5). As a consequence, there is an inertial form 
of (6.1), which reproduces all the dynamics of the Navier-Stokes equations (6.1). 
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6.2. The Barotropic Equations of the Atmosphere. 

All our previous results are also true after some modifications for the following 2D 
barotropic vorticity equation of the globe atmosphere, which was derived by Rossby in 
the 1930s, and was used by J.G. Charney, R. Fj0rtoft & J. von Neumann [CFN] in the 
first weather predictions in the 1950s: 

{ut + Vuit + 2fi cos 9k X u — uAu = /, 

where the space domain is S'^, the 2D sphere with radius a, A is the Laplace-Beltrami 
operator on 5"^, k is the unit outward normal vector, O is the angular velocity of the 
earth, and 9 is the colatitude of the earth. The existence and dimension estimates of 
attractors of these equations (in vorticity form) were studied in [W]. We remark here 
that all the results in Sections 1-5 hold true for equations (6.7). Particularly, we can 
obtain an inertial form of (6.7) with dimension given by an estimate similar to (5.21). 
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